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Abstract 



We consider the Harmonic crystal, a measure on R with Hamiltonian i2"(x) = 
J2i j ^ij'( x W — x (j)) + ^ Si( x W ~~ d(i)) 2 , where x, d are configurations, x(i), d(i) € 
M., i,j G Z d . The configuration d is given and considered as observations. The 
'couplings' Jj j are finite range. We use a version of the harness process to explicitly 
construct the unique infinite volume measure at finite temperature and to find the 
unique ground state configuration m corresponding to the Hamiltonian. 

Keywords Non homogeneous harmonic crystal, Harness process. 

1 Introduction 

The harnesses were introduced by Hammersley [10] to model the behavior of a crystal and 
to introduce a multi-dimension version of a martingale. Let P = (p(i,j), i,j G Z d ) be a 
homogeneous symmetric Markov transition matrix with p(i, i) = 0. A harness is a measure 
on M zd with the property that the mean height at i given the heights at all sites different 
of i is a P-weighted mean of the heights of the other sites. The serial harness is a Markov 
process on M z updated at all discrete times at all sites by the rule: substitute the height at 
site % by a P-weighted mean of the neighbors plus a centered independent random variable 
(the noise). Hsiao [TT] proposed a continuous-time version then called harness process in 
[7]. The heights are updated at Poisson epochs using the same rule as in the serial harness. 
If the noise is a centered Gaussian random variable, the reversible measure of the process 
is the harmonic crystal, that is, the Gibbs measure with Hamiltonian 

H(x):=Y,AMi)-x(3)) 2 (1-1) 

where x G K zd , x(i) for i G Z d represents the height at site i and Jij = p(i,j); the 
temperature is given by the variance of the noise. 

We study a version of the harness process with a external local data term. We can think 
that each site i G Z d has an additional "neighbor" with a fixed height d(i), the data. The 
updating of the height at i involves the data d(i) in the averaging. This is a "heat bath" 
dynamics associated to the quadratic (Gaussian) Hamiltonian (12.21) : at rate 1, the height 
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at site i is substituted by a random height distributed with the conditional distribution 
associated to the Hamiltonian, given the heights at the other sites. We show ergodicity 
of the harness process with data d G X, a set of configurations with a mild restriction 
on the growth, defined in (12. 5ft ; this extends the work of Hsiao [HI [12] who considered 
the case d(i) = 0. Ergodicity means that there exists a unique invariant measure for the 
process and that the process starting from measures concentrating mass on X converges 
to the invariant measure. The unique invariant measure is also reversible. We show that 
any infinite-volume Gibbs measure on X associated to the quadratic Hamiltonian (12.21) is 
invariant for the dynamics. This fact and ergodicity imply that there is only one infinite- 
volume Gibbs measure associated to the Hamiltonian. When the process is constructed in 
a finite subset A of Z d , the invariant measure is a harmonic crystal with external site-by- 
site field d. As A grows to Z d the harmonic crystals in A converge to the unique invariant 
measure for the harness process in Z d . If the data is flat (i.e. constant) then we are in the 
case of massive lattice models in quantum field theory which are well known (see [SI IS] ) - 

Shortly our method is as following. We slice the space-time configuration space S = 
M. zd x R + in pieces determined by the realizations of the Poisson processes governing the 
updating times. We show the convergence of the process almost surely for (almost) every 
slice. To this end we use Harris graphical construction of the process as a function of a 
space-time marked Poisson process of rate 1 on Z d x R. It is convenient to construct the 
process in an arbitrary interval [s,t] to be able to take the limit as s -> — oo which is 
equivalent to take the limit as t — > oo in distribution. At each Poisson epoch, the value of 
the process at corresponding site is substituted by an average of the values of the process 
at the neighboring sites plus an independent noise. When the construction is explored 
backwards in time the value of the process at site % at time t can be expressed in function 
of the probabilities of a random walk running backwards in time conditioned on the space- 
time epochs of the Poisson process. The walk is killed at a rate related to the weight of 
the external data and when it hits the boundary. The value of the process at site i at time 
t starting at time s is expressed as a sum of four terms: the contributions given by (1) the 
noise, (2) the external data, (3) the boundary condition (in case the process is studied in 
a finite region) and (4) the initial condition. The noise contribution is a martingale with 
uniformly bounded second moments so it converges as s — ► — oo. The data contribution 
converges to a deterministic function m, a harmonic function for a kernel associated to J 
and h. The boundary and initial contributions go to zero as the region grows to I* d and 
the initial time s goes to — oo, respectively. 

In Section [2] we state the main three theorems. The first one shows the existence of 
the Harness process in infinite volume. The second one says that the Harness process 
with external data is ergodic, that its unique invariant measure coincides with the infinite- 
volume Gibbs measure and that there is a unique infinite volume Gibbs measure for this 
Hamiltonian. The third theorem shows that the harmonic function m is the unique function 
minimizing the Hamiltonian with the data term (ground state). In Section [3] we construct 
the process and state intermediate results for the finite and infinite process. In Sections 
S] and we show that the finite process has as reversible measure the Gibbs measure in 
finite volume and that the space and time limits coincide. In Section [6] we show that the 
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harmonic function m is the unique minimizer. 

2 Main Results 

Let A C Z d be a finite set and define the Hamiltonian 

tf A (x) = Yl J ^( x « - x ^) 2 + ^E( x « " d ^)) 2 ( 2 - 2 ) 

x, d G M zd , Jjj- is a finite range pair potential, that is Jjj = if |z — j| > R > 0, d 
is a fixed configuration which can be taken as data and h > is a fixed parameter. In 
fact i^A = H\(-, d, /i) depends on d and h which are fixed and dropped from the notation 
unless necessary. Let xa be the configuration x restricted to the set A and the superposition 
configuration xaYa^ be given by x for sites in A and by y otherwise. Let |A| be the number 
of sites in A. The family of measures {/iA(-|yA c )' A C Z d , |A| < oo, y 6 l z }, defined by 

e -H A (x A y A c) 

AtA(dx A |yA=) = v /^(^a), (2.3) 

^A 

where /i A is the Lebesgue measure in M A and Z A is a normalizing constant, is called 
the specification associated to the Hamiltonian H. A measure /i defined on IR Z is said 
to satisfy the Dobrushin-Lanford-Ruelle (DLR) equations if its conditioned distributions 
coincide with the specifications: 

/i(dx A |x A c = y A c) = /i A (dx A |y A c) (2.4) 

for //-almost all y. 

We construct the dynamics in a set of configurations with limited growth. Let 

X = jx G R zd : \*(])W jl/R < oo for all i G Z d | (2.5) 

Recall that R is the radius of the interactions. We introduce the harness process {rjt(i), i G 
Z d } on X. Let a G [0,1], P = (p(i,j), i,j G Z d ) be a space homogeneous finite-range 
symmetric stochastic matrix (that is, p(i,j) > 0, ^2jP{i,j) = 1 for all i, p(i,i + j) = 
p(0,j), p(0,j) = if \j\ > R > and p(i,j) = p(j,i))- Symmetry is not necessary to 
define the process but it is natural in this context: we relate later p(i,j) with Jjj which is 
symmetric. 

The generator of the process acts on locally finite continuous functions / by 

L/(x) = Y WW = j °( dx ) If (« P * X + i 1 - «W d (*0 + *kx) - /(x)] , (2.6) 

keZ d k£Z d 
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where G(dx) = x2 / 2 dx is the standard Gaussian distribution (with zero mean and 

variance 1), called noise and P&x is the configuration defined by 

Wx)0-) = (^.f P(M)XW "> = *■ (2.7) 
I x(j) otherwise 

and efc(j) = l{/c = j}. In this process at rate 1 the height at site i is updated with 
a convex combination of the heights at the neighbors of i and the height d(i), plus an 
independent standard Gaussian variable. That is, if site i is updated at time t, the height 
at i is substituted by a random variable with the same law as 

a(P^)(i) + (l-a)d{i) + Z (2.8) 

where Z is a standard Gaussian independent of r] t _. The real variable (I2.8p has law 
//{.;}(• |(^_)/nc), the distribution of the {i} coordinate given the values of j] t - in {i} c as 
defined in ( 12.31) . Our dynamics coincides with the so called "heat bath". 

The following three theorems are the aim of this work. The first one shows the existence 
of the process. 

Theorem 2.1. Assume d £ X. There exists a Markov process (r] t ) on X with generator L: 

]im-E[f( Vt+u ) - f( Vt )\%] = Lf(rj t ), (2.9) 

tt^O u 

where % is the sigma algebra generated by (r] s , s < t), the past of r] s up to time t. 

Theorem 12.11 is a consequence of a general existence result of Basis jl] ; we provide here 
an alternative construction. 

Let S(t) be the semigroup associated to L defined by (S(t)f)(rj) = E(f(r)t)\r)o = rf). It 
acts on measures by j fd(uS(t)) = j S"(t)/dz/. Our second result says that the harmonic 
crystal is the invariant measure for the harness process when p(0, 0) = and establishes 
time and space limits. 

Theorem 2.2. Assume d £ X and p(0,0) = 0. (i) The following time and space lim- 
its exist and are identical. For any initial measure v concentrating on X, any boundary 
conditions y £ X and any increasing sequence A /* Z d , 

lim uS(t) = lim ua := U- (2.10) 

(ii) fi is reversible for for all k £ 7L d , in particular it is reversible for L. (Hi) \i is 
the unique measure in X satisfying the DLR equations for the specifications $2.3\) with 
Jij = ap(i,j) and h = 1 — a. 

The theorem implies that the harness process is ergodic in X: there exists a unique 
invariant measure \x and the process starting in X converges to \i. It also implies that 
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for any boundary conditions in X the thermodynamic limit is unique (absence of phase 
transition). 

The time convergence was proven by Hsiao [H] and [32] in the case d(fc) = 0. The 
space convergence is contained in Spitzer [13] and Dobrushin [I], see also Caputo [3]. Our 
approach permits to construct simultaneously (coupling) realizations of the measures in 
all finite boxes and the infinite volume measure in such a way that the convergence is 
almost sure. Notice however that the uniqueness result is a consequence of the fact that 
there is only one infinite-volume invariant measure for the process and that any measure 
compatible with the DLR conditions is invariant. The space thermodynamical limit is not 
used in the proof of uniqueness. 

The condition p(0, 0) = is necessary to guarantee that the measure fi\ is invariant for 
the generator L^. 

In Statistical Mechanics it is natural to extend the specifications (12. 3p to a family of 
measures defined by 

/4(dx A |y Ac ) = e -/WA(xAyA») ^ (2 n) 

Z A 

for j3 > 0; (3 is called the inverse temperature. We consider in detail the case (3 = 1; the 
other cases reduce to this one using /3i?A( x AyA c ) d, h) = -£/a( x aVa c ' ^ '> where xP(i) = 
\Z~j3x(i), etc. 

When j3 = oo the randomness vanishes and /i^ is interpreted as the measure concen- 
trating mass on configurations xa minimizing H\(^x\y\c) for finite A. When A = Z d , we 
denote H(x) = i? Z d(x), an infinite sum only formally defined. In this case we need to give 
a sense to the word "minimizing" . If x differs from x on a finite set of sites A, then the 
infinite sums defining H(x) and H(x) differ only on a finite number of summands. We 
define H(5i) — H(x.) as the difference of the corresponding different summands, that is, 
if A ( x A x A c ) — -^a( x a x a c )- We say that x minimizes H(-) if -ff(x) —H(~x) > for all x local 
modification of x. Measures defined on X compatible with the specifications (12.31) with 
j3 = oo are called ground states; see appendix B of van Enter, Fernandez and Sokal [6] for 
details. The ground states are concentrated on minimizing configurations. 

Let K(i,j) be the probability that the walk with rates aP (a walk killed at rate 1 — a) 
is killed at site j when starting at site i. In the next theorem we show that a delta measure 
concentrating mass in the configuration given by the K-average of d is the unique ground 
state. 

Theorem 2.3. Assume p(0, 0) = and that d £ X and let m be the configuration given 
by 

m(i) := 5^K(i,j')d(i) for all i £ Z d (2.12) 

3 

then m minimizes the Hamiltonian H and the delta measure concentrating mass on m is 
the unique ground state for specifications A2.3\) in X. 
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Using the Kolmogorov Backwards equation for the walk with transitions aP killed at 
rate 1 — a and boundary conditions d, we see that m satisfies the equation 

m(i) = <xp{hj)**U) + (1 - «) d «- (2-13) 

That is, m is a harmonic function for a transition matrix associated to P and a in an 
extended graph with "boundary conditions" d. The extended graph has vertices in 7L d U 
(Z d )* where (Z d )* is a copy of Z d , and edges G Z d xZ d : = l}U{(z',i*) : i G Z d 

and z* is the copy of i in (Z d )*}. The boundary conditions are fixed in (Z d )* equal to d; 
d(i) is the value of the boundary condition at i* and (1 — a) is the weight of the edge 
ap(i,j) is the weight of the edge 

3 Harris graphical construction 

The proof of the above theorems are based on an adaptation of the Harris graphical con- 
struction for the harness process proposed by the first two authors in [7] . The core of the 
construction is a rate-one space-time Poisson process KonZ d x R. This can be thought of 
as a product of homogeneous one-dimensional Poisson process in M, one for each i G Z d . 
Space-time points in N are denoted (i, r) and called epochs. To each (i, r) G 2NT we attach 
two independent marks: t) and ip{i,r), where £(«,t) is a Gaussian random variable 
with zero mean and variance 1 called noise and <p(i,r) are variables whose distribution is 
described later; these two families are iid and mutually independent and independent of 
N. We denote P and E the probability and expectation induced by these Poisson processes 
with marks. 

The Harness process is realized as a function of the Poisson epochs and the marks 
£ (for the moment we do not use the marks </?) as follows. The height at each site % 
only changes at times (i, s) G N. Assuming that the configuration at time s— is r] s _ and 
(i, s) G N is a Poisson epoch, then the configuration at time s at sites k ^ i does not 
change (r) s (k) = r) s -(k)) and 

7i 8 (i) = a{Pv a -)(i) + {l-a)d(i)+Z{i,s) (3.14) 

where Pi is defined in (12.71) . In other words, at the Poisson epoch (i,s), the height at % 
is substituted by a average of the other sites and the external data at i plus a Normal 
random variable independent of "everything" else. By a standard percolation argument 
this construction can be performed for small time intervals, so that Z d is partitioned in non 
interacting pieces); we sketch it in the proof of 13.41 in Section ?? later. We give another 
construction based on a "dual" representation using the variables tp. 
Let the variables tp(i, r) be independent with law 

P(tp(i,r) = j) = ap(i,j), (3.15) 
P(^(i,r)=0 = l-a, 
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where i* (G - Z d ) is the copy of % in (Z d )*. The parameter a G (0, 1) is later chosen as in 
Theorem 12.21 Define a family of backward random walks indexed by (7, t), the space-time 
starting point, as a deterministic function of the Poisson epochs X and the marks <p (here 
we do not use the marks £). Fix t G R and 7 G Z d U (Z d )*. For each s < i we define 
c?^ G Z d U (Z d )* as the position of a random walk going backwards in time with initial 
position (at time t) aj tt ^ =7 and evolving with the following rules. The walk does not 
move between Poisson epochs and for s < t, if at time s+ the walk is at site i G Z d and 
(i, s) G X, then at time s the walk jumps to the position tp(i, s): for s < t, 

7 jV(«,s) if 1] = % G Z d and (i, s) G X 



if ^ )t] = 1* G (Z d )* or aj +it] = i G Z d and (z, s) £ X 



We say that the walk a? t , is absorbed at j G Z d at time s if aj u t , = j* for u < s and 
aJ uA G Z d for u G (s,t\. The family ((0^, u<t), t G M,7 G Z d U (Z d )*)) is a function of 
X and but we drop this dependence in the notation. 

A key object in this analysis is the law of the walk a conditioned on a realization of 
the Poisson epochs X. For % G Z d let 

6 M (i,j) = PH,t]=i|X) (3.17) 

be the probability the walk a starting at % at time t to be at j at time u given the Poisson 
epochs. The probabilities fer Ujt i (i, j) are function of the Poisson epochs X and do noi depend 
on (p. Analogously, for (j, r) G IN" define the probability of absorption at time r at site j 
given the Poisson epochs by 

a [T ,t]^f) = P (4,*] = e Z "> for M e MM ( 3 - 18 ) 

Notice that 

a[r,t](hf) = (1 - a) 6[ T)t ](z, j) (3.19) 

Define the X-valued process 77^] in the time interval [s, t] with s < t and initial condition 
z at time s by 

= X] b lT,t](hMtii r ) + a ht](«,J*)d(j) + 6 [S)t ](i, J>(j), (3.20) 

(j,r)eMls,t] jGZd 

where X[s, i] = {(j, r) G X : r G [s, t]} and recall r) is the noise associated to the 
epoch (j, r). Under this construction 77ut](z) is a function of X[s, t] and the corresponding 
noises £; it is an average determined by X of the noises £, the external field d and the 
initial condition z at time s. Our goal is to prove the following 

3.1. For each set and z G X the process {rj[ s ,t], t > s) defined in \3.20) is "well- defined. 
Namely, the sums in $3. 2(A) are finite with probability 1 and t)\ s a G X for all s < t. 
Furthermore the process is Markovian with generator L given in ( QS| ) and initial condition 
z at time s. 
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3.2. For any configuration z G X and fixed t 6 1, the limit 

lim r][ s>t] (i) := rj t (i) (3.21) 

s — y — oo 

exists with probability one and does not depend on z. The process (r] t , t G M) is a stationary 
Markov process with generator L given in Ii2.6]) . 

3.3. Call fi the marginal law ofrjt (which does not depend on t), then \x satisfies the DLR 
equations (12.4| . 

3.4. If jl satisfies the DLR equations then fx is reversible for the process (fi[ s ,t], t > s). 

3.5. fi, the marginal law of rj t , is the unique infinite volume Gibbs measure on X for the 
specifications ( 12.31) . 

In other words, the strategy of our proof is to construct a stationary process in infi- 
nite volume whose time marginal is the unique Gibbs measure for the Hamiltonian (12.21) . 
Theorem 12.11 follows from 13.11 



Finite volume We start considering the Hammersley processes in finite volume A C 
Z d with boundary conditions y G X. The updates occurs at space-time Poisson epochs 
(i,s) G N as follows: 



(5>(*,J>s4j) + X>(i,.7')y0')) + (1 - a)d(i) +£M (3.22) 



j'gA jeA c 



That is, at the space-time Poisson epochs, the process substitutes the value at i by an 
average of the values at the other sites and the external data, the values outside A are kept 
fixed and given by the boundary configuration y. 

The construction of the finite harness process in A with boundary configuration y goes 
along the same lines as in infinite volume, the difference is that the probabilities b and a 
are computed for walks that are also absorbed at A c . 

Call A* the copy of A in (Z d )*. We define a family of backward random walks absorbed 
at A c U A* indexed by (7, [s,t], A), with space-time starting point (7, t) with 7 G A. The 
rules now are 

7 (<p(i, s) if = z e A and (z, s) G X 

% ' t] ' A H + ,t ]A if tfs + , t] ,A e (Z d )* U A c or ^ s+)t])A = i G Z d and (z, 3) i N 

The only difference is that now the walk is absorbed at A* and at A c . If the walk starts at 
A, it will be absorbed either at A* or at A c . Calling 

tt([M],A) = {(j,r) G X : r G [s,t),j G A}, 
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the family ((o7 ut -\ A , u G [s,t]), t G R, 7 G A is a function of N([s,£], A) and the associated 
tp, but we drop this dependence in the notation. For i 6 A and u < t define 

W(«'.i) = p «*]^ = jM. ( 3 - 24 ) 

the transition probabilities for the walk absorbed at and A c given the Poisson epochs. 
The probabilities of absorption at time r at site 7 G A* U A c given the Poisson epochs are 
defined by 

a[r,t]Ah 7) = p (4,t],A = 7, 4^],a G A ' for u e ( r ' *] M ( 3 - 25 ) 

We define the R A -valued process T)[ a ,t],A m the time interval [s,t] with s < t, initial 
condition z G X and boundary conditions y G X at time s by 

^[a,a],A = ZA 

Vls,t}A i ) = (hr^A^MU^) + 0[r,t],A(«,J*)d(j) + o [T)t]iA (i, fc) y(fc)J 

(j',r)6>f([s,t],A) fceA c 

+ ^ 6 M,A(^i)zA(j)> foriGA (3.26) 

jeA 

where N([s, i],A) = {(j,r) G IN" : j G A, r G [s, t]}, the value £(j, r) is the Gaussian 
random variable associated to the Poisson epoch (j, r) G X. 

With the above construction nr^iA^) is an average determined by N (through the 
weights a and b) of the noise £, the external field d, the boundary configuration y and the 
initial condition z at time s. We drop these dependences in the notation. 

We prove the following facts about the finite process: 

3.6. For each s G R the process (r)[ St t],A, t > s) is Markov with initial condition z A; boundary 
conditions va c and generator 

£a/(x a ) = Y L kf "(xaYaO 

fceA 

/• (3-27) 
= G(dx)[f([aP k (x A y A c) + (1 - a)e fc d(A;) + e fc x] A ) - /(x A )] , 

fceA ^ 

where f is a bounded function depending only on coordinates in A 

3.7. TTie measure /iA('|yA c ) grwen m f)2.3p is reversible for the process (^utiA) t > s) for 
each fixed s G R. 

3.8. For any configuration z G X and /zxed £ G R, tae Zzmit 

lim r][s,t].A(i) ■= »7t,A(i) (3.28) 

exists with probability one and does not depend on z. Tae process (77* a > i 6 R) « stationary 
with time-marginal /iA( - |yA c ) ? which is the unique invariant measure for the process. 
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3.9. For any configuration y G X and fixed t el, the limit 

lim VtA (i) = r) t (i) (3.29) 

A^Z d 

holds with probability one and does not depend on y. 
Theorem [275] follows from 13721 to 13791 . 

4 Finite volume 

In this section we prove [3761 13.71 and 13.81 for the finite process 77^. 

Generator Proof of \3.6l From (13.261) it follows that 
Vls,u+h]Ai) = l{|>T([«,« + /i],{i})|=0}?7 M ^(z) + l{|N([«,« + /i], {0)1 = 1} (4-30) 

x (oi ^2 r p(i, j)V[s, u ]Aj) + a ^ K*» J)y(i) + C 1 - a)d(z) + £(r, i)J + other terms 

jeA jeA c 

where r G (w, u + /i), {|N([rt, it + h], {i})\ = k} is the event "there are exactly k Poisson 
epochs in [u, u + h] x {2}" and the "other terms" are related to the presence of more than 
one Poisson epoch in U {%) x [u, u + h] which has probability of order h 2 , where U(i) is the 
cube centered at i with side R. The independence properties of the Poisson process and 
(I4.30p show that the process is Markovian with generator (I3.27p . □ 



Invariant measure Proof of \3. 7[ Let p(0, 0) = 0, < h < 1 and = (1 — h)p(i,j) 
and h = 1 — a. We prove that /ia defined by ( 12731 on M A is reversible for each one of the 
generators Lfc defined in (I3.27p . that is, for /, g : IR A — > R, 

Mfl'W) = HA{fL k g). (4.31) 

For a fixed configuration z and site fc let fc G A, let 

z(fc) = ^ J fc)i z(i) + hd(fc) (4.32) 

Let z = XAyAc A direct calculation yields 

Rk := Jk M^) ~ <k)f + h(z(k) - d(k)) 2 (4.33) 

i: i^k 

= (z(k) - m? + E J *.*(»(0 - z(*0) 2 + - d(fc)) 2 . (4.34) 

i: ij^k 

It follows that the conditional law at coordinate k given the heigths at the other coordinates 
is a Gaussian with mean z(k): 

m (dz|z { , }c ) = -L e -(-*( fe )) 2 dz (4.35) 
V 2ir 
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so that the updating at site k is done with the conditional distribution given the heights 
at the other sites. This implies reversibility; to show it we use the notations 

£ fe (z)=expj- JiMQ~*U)) 2 - E^W-dW) 2 ] 

E k (z, z{k)) = eJ-£ J k Mi) ~ <k)f - h(z(k) - d(k)) 2 I 

I i:i^k ) 

We represent the left-hand side in (14.311) as 

HA(gL k f) = sl - s 2 k , 

where 

Sl--=^r fl[dz(i)E k (z)E k (z,z(k)) f dxe~ x2 g(z)f ((1 - h)P k z + hd{k) + e k x) 

A J - gA JR. 

% --=Y~J \{^)E k ^)E k Mk))g{z)f{z) 
Using (14.341) . we obtain 

S * = Y~J II dz ( i )J5*(a)J5*(».S(*))«p{-(z(fc) - z(fc)) 2 } x 

dxe-'Vz)/ ((1 - /i)P fc z + /id(fc) + e fe x) 



Next we change the variables as follows 

v(*) = 

and 



z(/c) + x, if z = A;; 
v(i) = z(i), otherwise; 



y = z(k) - z(k). 

Observe that v(/c) = z(k). Therefore z(i) — z(k) = v(i) — v(k) for any i ^ k. Using 
first the above substitution and then relation ( 14. 34ft again yields 

S * = T /n dv (^ fc ( v )^(v,v(fc))exp{-(v(A;)-v(A;)) 2 } 

A ^ ieA 

x / dye- y2 g((v(k) + y)v A \ k )f(v) 
Jr 

= JT I T\^(i)E k {v)E k {v,v(k)) I dye-y 2 g((l-h)P k v + hd(k) + e k y)f(v). 
a J . gA Jm 

We thus obtain fi A (fL k g) = Sl - Sl = fi A (gL k f). □ 
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Limiting Distributions Proof of \3.S\ We consider separately the four terms in (13.261) . 

Define 

A M ,a«:= E hrAAhjWiT) (4.36) 
0>)eN([ Sl t],A) 

B [s , t] , A (i) := a \r,t]AhfW) (4-37) 
0-,r)eN([ s ,t],A) 

c m ,a« : = E E a M,A(*> fc ) ( 4 - 38 ) 

(j»eW([s,i],A) fceA c 

^M,a« :=^&M],a(m>(j) (4-39) 

Since the variables £ are independent of the variables b, (A\t- u j],A(i) , u > 0) is a martingale 
for the sigma field 5F U generated by N([t — w, £], A) and the associated Gaussian variables £. 
The variance at time u is given by 

E(A [t _ uAA (i)) 2 = E E (&M],A(i,i)) 2 (4.40) 

(i,T)6N([i-«,t],A) 

by conditioning on the Poisson epochs and deleting the cross terms by independence of the 
different £. The last expression is bounded by 

< E E = E f e- {1 - a)r Ql^3)dr < ^— (4.41) 

uniformly in A, where Q r A (i,j) is the probability for a continuous time random walk with 
rates p absorbed at A c starting at i to be at j by time r. The factor e~^~ a ^ r corresponds 
to the probability that the walk is not killed in the time interval [t, t — r]. The Martingale 
convergence theorem implies that 

lim A [t _ u>t]A (i) = A tA {i) a.s. (4.42) 

u— >oo 

The process (A tA , t e M) is Markov, stationary with generator 

E ^/(xaYac) = E / G ( dx ) [/(«^'( x Ay^ + e k x)) - /(x A y Ao )] , (4.43) 

and E(A tA (i)) 2 <(l-ay\ 

We show now that the limit in (I4.37P is given by 

lim B [t _ uAA (i) = VK A (i,i)d(i) := m A (z) (4.44) 

ieA 
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where K\(i,j) is the probability the random walk to be killed at site j, if j G A, or to be 
absorbed, if j G A c . As in (12.131) . niA satisfies 

m A (i) = y~] ap(i, k)m A (k) + (1 - a)d(i). (4.45) 

fceA 

For v > u, the following "Markov property" holds 

B[t- v ,t],A{i) '■= -B[t-u,t],A(«) + / b[t-u,t],A (z, k)B [t _ Vtt _ u]A (k) (4.46) 

fceA 

Applying (14.451) to each Poisson epoch (j, r) we have 

m A (z) = £[t-«,t], A (i) + 5^ 6 [t-u,t].A(i» fc)m A (fc) (4.47) 

fceA 

which implies that m A is invariant for this dynamics. For fixed s the process (B\ St t],A> t > s) 
on M A is Markov with generator 

£>£/(x A ) = J2 [/([«Pfc(x A o Ac ) + (l-«)e fc d(fc)] A ) -/(x A )], (4.48) 

fceA fceA 

where the o is the "all zero" configuration: o(i) = 0. Finally, using (14.471) . 
|£[i- w ,t],A(X) -m A (z)| < 

{i,k)\m A (k)\ (4.49) 

fc 

By (14.411) the right hand side converges to zero as u — > oo uniformly in A, proving (14.441) . 
A similar argument shows that 



lim C [t _ U)t]>A (i) = V K A (z, j)y(j) := r A (i) (4.50) 

u— >oo » 

jeA c 

Finally, by pi) . 

lim IVm,aW = (4.51) 

u— +oo 

The limits (1Q2]h (jOi) . fT4T50l and flCTj) show (ET2S1) . The resulting limit is 

r)t,A = A tA + m A + r A (4.52) 

By construction the law of this limit does not depend on t. One proves that the process 
{Vt,Ai ^ G M) is Markov like in 13.61 Since the limit (I4.52p does not depend on the initial 
configuration z, this shows that the process (^,a) has a unique invariant measure. Since 
by 13.81 /i A (-|y A c) is invariant for this process, for each t e 1 the marginal law of r]t t \ is 

A*A(-|yA<0- □ 
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5 Infinite volume and thermodynamic limit 



Existence of the process Proof of \3.1\ The fact that the sums in ( 13.201) are finite 
follows immediately from the finite range condition on p. The proof that the dynamics is 
the harness process follows then from an argument similar to the one in 13.61 □ 

Limiting stationary processes Here we prove 13.21 and 13.91 To prove 13.21 we consider 

(5.53) 



in (|3.20|) 


Define 


A M(i) ■ 


Y b lT,t](hMUi T 




0» g :n-(M) 


B [;t]{i) ■ 


Y a [T,t]{i,j*)d(j) 











(5.54) 
(5.55) 

Lemma 5.1. Under the hypothesis of Theorem \2.2l 

lim A [t _ U}t] (i) = A t (i) (5.56) 

u— >oo 

and (A t , t G R) is a stationary Markov process with generator 

YL£f(*) = Y f G(dx)[f(aP k x + e k x)-f(x)], (5.57) 

k k J 

with E(A t (i)) 2 < (1 — a)" 1 . Furthermore, 

lim A tA {i) = A t (i) (5.58) 

KfZ d 

Proof. As in the proof of (14. 42 p . (A[ t _ u t ](i), u > 0) is a martingale with uniformly bounded 
second moments. Hence the limit (15.560 exists. Since b[ T j],A(i, j) S b[ T ,t](h j), 

\A t (i) - A t>A (i)\ < Y ( b [r,t](hj) - & M,a(«, j)M{j,r)\ (5.59) 

(j,t)6X1-oo,*] 

The sum is finite because the variance of A t (i) is bounded and it converges to zero by 
monotone convergence. 

The process (A t , t G M.) is stationary by construction. It is Markov with generator 
(I5.57P as a consequence of (I5.58P and 13.61 □ 

Lemma 5.2. //d G X then |m(z)| < oo for all i G Z d and 

lim B [t _ ut ](i) = m(i) (5.60) 

u— >oo 

and 

lim mA(i) = m(i) (5.61) 

A/1 d 
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r 

Proof. Remark that K(i,j) < al * J because at least 



R 

achieve the point j from i and during the time erf is not absorbed. 



1 jumps are needed to 



oo,t] 

m(i)| <^2K(i,j)\d(j)\ < $>[^]|d(j)| < oo (5.62) 



by the definition of X. The proof of (15.601) works as the proof of fl4.44[) . 
On the other hand K\(i,j) f K(i,j), hence 

|m(z) - m A (i)| < Y^{K(i,j) - K A (i, j))\d(j)\ - (5.63) 

j 

by monotonic convergence. □ 
Lemma 5.3. Ify G X then for r A defined in (14.501) . 

lim r A (i) = (5.64) 
Proof. As in the proof of the previous lemma, 

Mi) | < ^«^|y(j)l (5-65) 

jeA c 

which converges to zero as A /* Z d because the sum is finite by definition of X. □ 

Proof of \3.S\ Calling rj t =: A t + m, the convergence follows from (I5.56p . (I5.60p and (14.51 j) 
(which holds uniformly in A). □ 

Proof o/rO Since 

Vt,A = A,A + m A + r A and r) t = A t + m , (5.66) 
the convergence follows from ( 15. 58ft , 15.611 and 15.641 □ 

Limiting process satisfies DLR conditions Proof of \3.3\ . Here we use the thermody- 
namic limit. The measures yU A (-|y A c) are compatible with the DLR conditions for subsets 
of A. By (I3.29P the measure /i is the limit as A increases to Z d of /iA("|yA c ), independently 
of y G X. Hence, also \x is compatible with the DLR conditions. □ 
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DLR measures are invariant Proof oj )3.4\ This proof does not use the thermodynamic 



limit. We first prove that Jl is invariant. Let rj be a configuration sampled with /L Let 
r}(k) be defined in (I4.32p . Let Z be a standard Gaussian variable independent of 77. Since 
Jl satisfies the DLR equations, a computation like in (I4.34p shows that the variable 

r](k) + Z has law /5(x(A;) G ■ | X{fc}c = rj^y) 

and the configuration 

9(k, t], Z) = t] - 5 v{k) + $rj( k )+z 

obtained by substituting the value at k by r} k + Z has law Jl. The Harness dynamics does 
the same substitution with the Gaussian noises £ at the updating epochs; see (I3.14p . This 
means that after each Poisson epoch, the distribution of the updated configuration is the 
same as the one just before the updating. But since there are infinitely many sites, there is 
no "first Poisson epoch" to apply the rule and show 13.41 directly. To overcome the difficulty 
consider a time t small enough such that 7L d is partitioned in finite sets A^, t = 1,2, .. . 
(depending on the Poisson epochs in the interval [0,t]) in such a way that if there is an 
epoch at time s G [0,t] at site i G Ag, then {j G Z d : \\i — j\\ < R} C A^, where we recall 
R is such that p(0,j) = if ||j|| > R. A standard percolation argument shows that this is 
possible. 

Enumerate the epochs as follows. Start ordering cronologically the epochs in Ai; con- 
tinue with the epochs in A 2 and so on. Under this ordering, if n(i, s) is the label of epoch 
(z, s) then n(i, s) < n(i', s') if s < s' or if i G A^, i' = A^ and £ < £'. Call (i n , s n ) the nth 
Poisson epoch in this labeling and £ n = £(z n , s n ) the associated Gaussian variable. Choose 
77 with law Jl and define inductively t] = rj and for n > 1, 

By the previous considerations r] n has law Jl for all n. By the definition of A^, the up- 
dating of sites in A^ in the time interval [0, t] depends on the initial configuration rj only 
through the values at sites in A^. On the other hand, for any cylinder function with sup- 
port on a finite set A, there exists an n(A, 3\f) < 00 such that A C ■ ■ ■ ,i n (A,jf)} an d 
j Jl(dr])E(f (ri n (A ^)) I N) = Jlf for almost all 3\f. Hence 

JlS(t)f = j /I(d7 7 )E(E(/(7 7n(A ^)|K)) = E(| Jlidr)) E(/(77 n(A>w) ) I >[)) = Jif 

This shows invariance of Jl. 

A similar argument shows reversibility. Going backwards in time, the law of r] n _i(i n ) — 
fj n {in) is a standard Gaussian variable independent of r] n . The construction can be done 
using the same Poisson epochs and partition (A^). □ 



6 Zero temperature 

Proof of Theorem We need to show that the configuration m defined in (I2.12p is com- 
patible with the specifications (12.31) at zero temperature. Let A C Z d be a finite volume. 
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We want to show that rriA given by f)2.12p is the xa which minimizes 

#(x A m A c) = a E p(m)( x (0 - x (i)) 2 + « E /c )( x ( i ) ~ m (^))' 



(6.67) 



ijgA 



i€A 

fceA c 



+(1 - a) J>(i) - d(i)) S 



ieA 



The Hamiltonian (I6.67P can be rewritten as 



iJ(x A m Ac ) = «^^^,A;)(x(2)-m(A ; )) 2 + (l-«)^(x( ? )-d(0)^ 

«eA fcgzrf 



ieA 



+ 2a E P{iJ) [ ~ ( X W ~ m 0'))( x 0') - m (0) 



(6.68) 



+ (x(i) - m(j))(m(j) - m(z')) - (x(j) - m(z'))(m(j) - m(z')) + (m(j) - m(z')) : 



then 



Let rj = m(z) 2 - |a J]fcez d ^) m (^) 2 + (1 - a)d(i) 2 

i7(x A m A c) = 5^(x(z) - m(z')) 2 - E r * 

ieA ieA 

+ 2a E p(z, j) [ - (x(i) - m(j))(x(j) - m(z)) + (x(i) - m(j))(m(j) - m(i)) 

- (x(j) - m(z))(m(j) - m(i)) + (m(j) - m(z)) 2 



ijeA 



To find the configuration x A minimizing (I6.67P we take the derivatives for ieA: 
dH 



<9x(z) 



2(x(i) - m(i)) 
+ a E p(z, fc) - 2(x(A) - m(z)) + 2(m(A;) - m(z')) 



fceA 



Then we obtain the following system of the equations 

x(z) — a ^^p(z, fc)x(fc) = m(z) — a ^^p(z, A;)m(/c). 



(6.69) 



fceA 



fceA 



It is clear that x(z) = m(z') is a solution of (I6.69p . Since if A is convex, m is the minimum. 
Since if A is a second order polynomial the minimum is unique. This shows that m belongs 
to the support of a ground state measure. 

Uniqueness Let m 6 X be a minimizer of H. Then fn(z') = a ^jP(z, j)m(j) + (1 — 
a)d(j) minimizes if(x{j}m{j} C ) for each i G Z d . So m is invariant for the dynamics 



E L ^( x ) = E [f( aP ^) + (1 - «)e*d(*)) - /(x) 

fcez d fcez d 



(6.70) 
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because this dynamics chooses a site i at random times and substitutes the height at i with 
the value minimizing H ('x^m^yc) . The graphical construction of this process with initial 
measure m and the invariance of m gives 

m(i) = B [t _ u fi{i) + ^ h-u,t](h k)m(k) (6.71) 

kez d 

Subtracting this equation from (I4.47P with A = Z d we get 

|m(i) - m(z') | < h-u,t] (i, fc) |m(fc) - m(A;) | (6.72) 

fcez d 

This is summable because both m ad m belong to X and since lim^oo b\t- u ,t] {i, k) = 0, 
then m = m. □ 
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